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Abstract 

The  transient  behaviour  of  a  solid  oxide  fuel  cell  (SOFC)  caused  by  load  changes  is  investigated.  For  this  purpose  a  step  function  of  the 
current  density  is  set  to  the  fuel  cell  and  the  response  of  the  cell  voltage  is  calculated  by  applying  a  special  SOFC  simulation  code.  The 
relaxation  time  is  strongly  related  to  the  transient  temperature  distribution  of  the  solid  cell  structure.  Therefore,  it  is  in  the  order  of  some 
minutes  depending  on  the  design  parameters  and  the  operating  conditions.  As  a  remarkable  result  the  cell  voltage  undershoots  beyond  the 
new  steady-state  value  from  where  the  relaxation  process  starts.  This  effect  is  discussed  and  explained.  The  transient  process  is  generalized 
by  introducing  two  dimensionless  groups,  the  modified  Fourier  number,  Fo,  and  the  source  term  number,  So.  In  terms  of  these  characteristic 
numbers  a  simple  relationship  is  derived  to  determine  the  relaxation  time  of  arbitrary  cell  designs  under  arbitrary  operating  conditions. 
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1.  Introduction 


Load  change  is  a  frequently  occurring  transient  process 
during  solid  oxide  fuel  cell  (SOFC)  plant  operation.  The 
control  system  has  to  react  on  such  event.  Therefore,  knowl¬ 
edge  about  the  response  time  of  the  stack  is  required. 

A  simulation  code  as  described  in  Refs.  [1,2]  was  applied 
to  calculate  the  transient  behaviour  of  the  SOFC  stack.  It 
predicts  the  operating  conditions  for  various  design  parame¬ 
ters  in  a  spatial  and  temporal  resolution. 

The  time-dependent  change  of  the  cell  voltage  caused  by 
a  jump  function  set  to  the  current  density  is  used  as  a  criterion 
for  the  determination  of  the  relaxation  time,  A  t.  It  is  defined 
as  the  period  necessary  to  recover  90%  of  the  dynamic  voltage 
drop  resulting  from  the  difference  )  -  U(t-  +0). 

For  this  study  numerous  transient  computations  have  been 
performed.  In  order  to  cover  a  broad  range  of  future  appli¬ 
cations  the  design  and  operating  parameters  have  been  varied 
widely.  In  particular,  the  dimension  of  the  bipolar  plate,  the 
channel  height  and  the  properties  of  the  solid  material  were 
changed  by  one  order  of  magnitude.  Furthermore,  various 
assumptions  for  the  fuel  utilization,  the  stoichiometry  factor 
and  the  current  density  have  been  applied.  Finally,  calcula¬ 
tions  without  and  with  internal  reforming  have  been  con¬ 
ducted.  In  the  latter  case  the  fuel  was  a  30%  pre-reformed 
gas  resulting  from  a  H20/CH4  =  2.5  mixture  at  CO  shift- 
equilibrium.  The  corresponding  gas  composition  is: 

0378-7753/95/S09.50  ©  1995  Elsevier  Science  S.A.  All  rights  reserved 
SSDI 03 7 8 -775  3  (95  >02263 -4 


CH4=  17.10%,  H2  =  26.26%,  CO  =  2.94%,  C02  =  4.36%, 
H20  =  49.34%. 

The  only  parameters  to  be  unchanged  are  the  active  elec¬ 
trochemical  area  of  100  mm  X  100  mm,  the  number  of  chan¬ 
nels,  ic  =  20,  and  the  rib  width  bT  =  2  mm.  In  Fig.  1  a  sketch 
of  the  cell  design  is  depicted. 

The  computation  of  transient  processes  is  CPU-time  con¬ 
suming,  and  not  all  design  engineers  have  an  appropriate  code 
available.  Therefore,  an  attempt  was  made  to  facilitate  the 
prediction  of  the  relaxation  time  by  developing  a  simple  algo¬ 
rithm.  This  is  the  subject  of  the  second  part  of  this  contribu¬ 
tion. 


2.  Computational  results 

In  order  to  emphasize  transient  effects  higher  steps  of  the 
current  density  are  applied  in  this  study  than  usually  will 


E.  Achenbach  /  Journal  of  Power  Sources  57  (1995)  105-109 


106 


Time  [.] 


Fig.  2.  Transient  cell  voltage  during  load  change  from  j,  to  j2  =  500  mA/ 
cm2;  ceramic  bipolar  plate. 
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Fig.  3.  Transient  cell  voltage  during  load  change  from  j,  to  j2=  500  mA/ 
cm2;  metallic  bipolar  plate. 


occur  during  normal  fuel  cell  operation.  For  this  reason  the 
cell  must  be  supplied  in  advance  with  that  amount  of  fuel 
which  is  required  to  achieve  the  new  high  performance  at  a 
given  fuel  utilization.  Those  calculations  simulating  a 
decreasing  current-density  step  start  with  a  given  high  utili¬ 
zation  and  reach  a  lower  value  for  the  new  steady  state. 

Figs.  2  and  3  show  the  evolution  of  the  cell  voltage  as  a 
function  of  time  for  cells  with  ceramic  or  metallic  bipolar 
plates,  respectively.  The  configuration  of  both  designs  is  iden¬ 
tical.  For  t  <  0  the  cell  voltage  corresponds  to  the  initial  cur¬ 
rent  density,  7).  The  asymptotic  value  for  t-*so  is  the  new 
steady-state  cell  voltage  for  the  current  density  j2.  The  inter¬ 
mediate  period  is  characterized  by  an  undershooting  of  the 
cell  voltage,  below  the  new  steady-state  value.  The  amount 
of  the  voltage  drop  correlates  with  the  magnitude  of  the  cur¬ 
rent-density  step.  The  relaxation  time,  however,  is  independ¬ 
ent  of  it. 

The  reason  that  the  voltage  indicates  the  undershoot  effect 
is,  that  immediately  after  the  jump  (t=  +0),  the  cell  tem¬ 
perature  is  still  low  and  related  to  the  initial  current  density. 
The  higher  amount  of  waste  heat  produced  at  j2  =  500  mA/ 
cm2  causes  the  cell  temperature  to  increase  which  results  in 
an  increase  in  the  cell  voltage  due  to  the  decreasing  internal 
cell  resistance. 

Figs.  2  and  3  demonstrate  that  the  relaxation  time  depends 
on  the  design  parameters  of  the  cell,  such  as  the  solid  material 


properties:  the  thermal  conductivity,  As,  the  heat  capacity,  cps, 
and  the  density,  ps.  Further  parameters  are  the  mass  and  the 
configuration  of  the  solid  structure  and  the  waste  heat  pro¬ 
duced.  The  ceramic  cell,  for  instance,  (As  =  2  W/(m  K), 
cps  =  400J/(kgK),ps  =  6600kg/m3)  needs  about  120  s,  the 
metallic  ( As  =  27,  cps  =  500,  ps  =  7800)  around  300  s  to  com¬ 
pensate  for  90%  of  the  dynamic  voltage  drop. 

Fig.  4  demonstrates  that  the  temporal  variation  of  the  cell 
voltage  is  closely  related  to  the  solid  material  temperature. 
The  figure  exhibits  the  maximum  and  minimum  transient 
temperatures  of  the  ceramic  cell.  The  temperatures  increase 
after  the  current-density  jump  and  approach  the  value  relevant 
for 72  =  500  mA/cm2.  The  response  time  of  7"max  nearly  cor¬ 
relates  with  that  of  the  cell  voltage  while  the  response  time 
of  Tm j,,  is  about  2.5  times  greater.  The  reason  for  this  is 
associated  with  the  counter-acting  effects  either  of  the  endo¬ 
thermic  steam  reforming  reaction  or  the  exothermic  electro¬ 
chemical  process  as  discussed  below. 

Fig.  5  shows  the  transient  cell  voltage  as  it  may  occur 
during  operation  of  a  large-scale  power  plant.  A  jump  of  the 
current  density  by  10%,  from  270  to  300  mA/cm2  is  assumed. 
Now  the  waste  heat  production  is  lower  as  compared  to  pre¬ 
ceding  cases  at  500  mA/ cm2.  This  results  in  a  larger  response 
time  of  Ar=220  s.  The  undershooting  voltage  drop  is  only 
6  mV. 
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Fig.  4.  Transient  minimum  and  maximum  solid  temperature  during  load 
change  from  j,  to  j2  =  500  mA/cm2;  ceramic  bipolar  plate. 
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Fig.  5.  Transient  cell  voltage  during  load  change  from  j,  =  270  to  j2  =  300 
mA/cm2.  Ceramic  bipolar  plate. 
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Fig.  6.  Transient  cell  voltage  during  load  change  from  j,  =  300  or  j,  =  500 
to  ;'2  =  400  mA/cm2. 

Finally  Fig.  6  exhibits  the  transient  cell  voltage  for  the 
situation  where  the  step  of  the  current  density  is  +  A j  and 
-A j,  respectively.  Once  the  current  density  is  set  from 
j1  =  300  mA/cm2,  once  from  yj  =  500  mA/cm2  to  the  new 
steady-state  value  of  y'2  =  400  mA/cm2.  Though  the  curves 
are  not  symmetric  with  respect  to  the  horizontal  line  through 
the  new  steady-state  voltage  value,  the  response  times  are  the 
same. 


3.  Evaluation  of  the  computational  results 

In  order  to  adjust  the  determination  of  the  relaxation  time 
to  the  design  engineer’s  requirements  an  attempt  was  made 
to  facilitate  the  prediction  procedure  by  introducing  a  simple 
algorithm.  For  this  purpose  the  governing  parameters  must 
be  identified. 

Since  all  chemical,  fluid  dynamic  and  mass-transfer  proc¬ 
esses  are  fast  compared  with  the  dissipation  of  heat  in  the 
stack,  the  dynamic  term  must  only  be  considered  in  the  energy 
balance  of  the  solid  structure.  This  means  that  a  quasi  steady- 
state  solution  of  the  problem  is  achieved  at  each  time  step 
except  for  the  temperature  of  the  solid  phase. 

Eq.  ( 1 )  represents  the  energy  balance  of  the  solid: 

4>rs  32rs  a  2ts 

(pCp)s_aT  _Aeff  *  ax2  +Aeff-V  <fy2 
a  2r 

+  4eff,z-^7  +&  (l) 

where  Ts  is  the  local  solid  temperature.  The  thermal  conduc¬ 
tivities,  Aeff,  of  the  quasi-homogeneous  stack  structure  are 
effective  quantities  which  must  separately  be  determined  for 
each  of  the  spatial  coordinates  according  to  Ref.  [  1  ] .  Since 
the  effective  thermal  conductivity  includes  the  contributions 
of  heat  conduction  through  the  solid  and  gaseous  phase  and 
of  radiation  it  is  a  non-isotropic  and  strongly  temperature- 
dependent  property. 

The  sum  of  the  heat  rates  produced  or  consumed  in  the 
chemical  and  electrochemical  processes  and  those  transferred 
from  the  fuel  and  the  air  to  the  walls  occur  as  the  source 
term,  <P. 


Eq.  ( 1 )  can  be  written  dimensionless  by  introducing  the 
following  quantities: 


where  A  T=TQ-TE  is  the  difference  of  the  outlet  and 
entrance  temperatures  of  the  cooling  air  stream.  The  quantity 
Ar  means  a  characteristic  time  interval  which  is  in  the  present 
case  the  relaxation  time.  As  a  characteristic  length  scale  the 
effective  thickness  of  the  solid  structure,  hcK,  occurs,  defined 
in  Eqs.  (3)  and  (4): 

hc„=(\-e)(dPEN+db  +  h^  +  h()  (3) 

dpEN(pcp)PEN-l-(pcp)b^h-l-  (*«+Af)J 

*  (pcp)PENdPEN  +  (pcp)b(db  +  /za-l-^f) 

(4) 

The  explanation  of  the  symbols  can  be  read  from  Fig.  7. 

The  total  specific  waste  heat,  Pw,  to  be  removed  from  the 
cell  can  be  expressed  by  the  electric  specific  output  P0  of  the 
cell  and  the  efficiency  tjc: 

Pw=Pc(i-Ve)/ve=Uj(l-Ve)/Vc  (5) 

Eq.  (1)  yields  together  with  Eqs.  (2)  and  (5)  the  dimen¬ 
sionless  form  of  the  energy  balance  of  the  solid: 

36>  AsAr  /Aeffxa2e>  1  — tfc  Uj  \ 

3t  (pcp) shlff  \  As  8£2  T)c  XAT/hJ 

The  solution  of  Eq.  (6)  depends  on  two  dimensionless 
groups.  The  first  is  the  Fourier  number,  Fo,  formed  with  the 
effective  height,  hc({,  of  the  cell  structure. 

Fo  =  AsAr/[(pcp)s/i2ff]  (7) 

The  second  group  derives  from  the  source  term  and  is  called, 
therefore,  the  source  term  number,  So: 

So  =  Uj  ( 1  -  t}e)  /  ( TjeAsA  T/hc{{)  ( 8 ) 

Its  physical  meaning  is  the  ratio  of  two  heat  fluxes,  i.e.  the 
waste  heat  flux  compared  with  a  heat  flux  resulting  from  heat 


Fig.  7.  Sketch  of  the  SOFC  solid  structure. 
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Fig.  8.  Fourier  number,  Fo,  vs.  source  term  number.  So.  Computational 
result  and  approximation  algorithm:  (  +  )  start-up  procedure. 
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Fig.  9.  Transient  cell  voltage  during  start-up  procedure  of  a  ceramic  cell. 
Current  density  step  from  j,  =  0  to  j2  =  300  mA/cm2. 

conduction  through  a  plate  when  /tcff>  is  its  thickness,  As  its 
thermal  conductivity  and  A  T  the  temperature  difference. 

As  mentioned  above  numerous  calculations  of  load  change 
processes  have  been  performed  to  reveal  the  relationship 
between  Fo  and  So.  For  this  purpose  the  relevant  parameters 
such  as  (pcp)s,  As,  hcfl,  AT  and  Uj  have  been  varied  within 


one  order  of  magnitude  to  cover  a  wide  range  of  applications. 
The  result  is  shown  in  Fig.  8.  The  data  can  be  approximated 
by  a  simple  equation: 

Fo  =  0J2Sol  1  (9) 

Since  the  So  number,  Eq.  (8),  only  contains  known  design 
data,  the  Fo  number  can  be  determined  from  Eq.  (9)  and 
hence  the  relaxation  time.  At. 

The  start-up  procedure  of  a  stack  appears  as  a  special  case 
of  the  load  change.  Here,  the  current  density  starts  from  =  0 
and  is  instantaneously  set  to  an  amount  of  j2  which  is  chosen 
in  the  present  example  to  be  j2  =  300  mA/cm2,  for  instance. 
Fig.  9  shows  the  temporal  development  of  the  cell  voltage 
for  a  ceramic  cell  heated  up  to  a  temperature  of  900  °C  before 
the  current  density  j2  is  demanded.  Within  about  200  s  the 
cell  voltage  reaches  the  maximum  followed  by  a  minor  grad¬ 
ual  decrease  until  the  final  steady  state  is  established  after 
about  900  s.  90%  of  the  voltage  difference  is  recovered  after 
a  time  interval  of  Ar=  150  s  which  is  the  relaxation  time  in 
the  present  definition.  This  result  expressed  in  terms  of  Fo 
and  So  numbers  is  indicated  in  Fig.  8  by  the  cross  symbol. 

The  development  of  the  whole  temperature  field  during  the 
start-up  procedure  is  given  in  Fig.  1 0.  The  particular  pictures 
for  the  period  up  to  160  s  indicate  a  continuous  change  of  the 
isotherm  pattern.  For  an  elapsed  time  greater  than  160  s  the 
change  is  no  longer  essential  and  a  uniform  temperature 
across  the  cell  is  nearly  established.  Only  the  maximum  and 
particularly  the  minimum  are  observed  to  be  more  pro¬ 
nounced  (see  also  Fig.  1 1)  at  longer  relaxation  times.  Due 
to  the  effect  of  the  slowly  relaxing  temperature  minimum  the 
average  temperature  decreases  slightly  resulting  in  a  minor 
diminution  of  the  cell  voltage. 
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Fig.  10.  Transient  temperature  field  during  start-up  procedure  of  a  ceramic  cell.  Current  density  step  from  j,=  0  to  j2  =  300  mA/cm2.  Fuel  from  left,  air  from 
below. 
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Fig.  1 1 .  Transient  minimum  and  maximum  solid  temperature  during  start¬ 
up  procedure  of  a  ceramic  cell.  Current  density  step  from  j,  =  0  to  j2  =  300 
mA/cm2. 


4.  Conclusions 


The  transient  behaviour  of  a  solid  oxide  fuel  cell  during 
load  change  is  treated  theoretically.  The  relaxation  time  is 
calculated  for  the  situation  where  an  immediate  step  function 
is  set  to  the  current  density.  The  response  of  the  cell  to  the 
increase  in  performance  is  discussed  with  a  view  to  the  tem¬ 
poral  development  of  the  cell  voltage.  The  computational  data 
show  that  an  undershoot  of  the  cell  voltage  occurs  followed 
by  a  relaxation  process.  Some  minutes  are  necessary  to  estab¬ 
lish  new  steady-state  conditions.  This  relaxation  process  is 
closely  related  to  the  transient  temperature  distribution  of  the 
solid  cell  structure,  since  the  internal  cell  resistances  are 
strongly  temperature  dependent.  Thus,  the  relaxation  time 
depends  on  the  thermal  properties,  size  and  configuration  of 
the  cell,  and  on  the  operating  conditions. 

An  attempt  is  made  to  facilitate  the  prediction  of  the  relax¬ 
ation  time  by  introducing  two  characteristic  numbers.  These 
are  the  modified  Fourier  number,  Fo,  and  the  source  term 
number,  So.  In  terms  of  these  dimensional  groups  the  follow¬ 
ing  simple  relationship  is  proposed  to  solve  this  problem. 

Fo  =  0.72So~" 

The  start-up  procedure  is  recognized  to  be  a  special  case 
of  the  load  change  process. 


5.  List  of  symbols 

b  width  (m) 

cp  heat  capacity  ( kJ/  ( kg  K) ) 

d  thickness  (m) 

Fo  Fourier  number,  Eq.  (7) 

h  height  (m) 

j  current  density  (A/ m2;  mA/ cm2) 

Pc  specific  electric  power  ( W  /  nr ) 

Pw  specific  waste  heat  flux  (W/m2) 

So  source  term  number,  Eq.  (8) 
t  time(s) 

T  temperature  ( K) 

U  cell  voltage  (V) 

x,  y,  z  coordinate  (m) 

Greek  symbols 

e  porosity 

©  dimensionless  temperature 

£,  17,  £  dimensionless  coordinates 
tjc  electric  efficiency  related  to  upper  heating  value  and 
fuel  converted 

A  thermal  conductivity  ( W/  ( m  K) ) 

<P  source  term  (W/m3) 

p  density  (kg/m3) 

r  dimensionless  time 

Subscripts 

b  bipolar  plate 

c  gas  channel 

eff  effective 

E  entrance 

O  outlet 

PEN  PEN  structure  (anode,  electrolyte,  cathode) 
s  solid  structure 
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